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and nonnegativeness constraints for solutions. 

Anton V. Proskurnikov, Member, IEEE , Ming Cao, Member, IEEE. 


Abstract —In the recent paper by Hamadeh et al. (2012) an elegant 
analytic criterion for incremental output feedback passivity (iOFP) of 
cyclic feedback systems (CFS) has been reported, assuming that the 
constituent subsystems are incrementally output strictly passive (iOSP). 
This criterion was used to prove that a network of identical CFS can 
be synchronized under sufficiently strong linear diffusive coupling. A 
very important class of CFS consists of biological oscillators, named 
after Brian Goodwin and describing self-regulated chains of enzymatic 
reactions, where the product of each reaction catalyzes the next reaction, 
while the last product inhibits the first reaction in the chain. Goodwin’s 
oscillators are used, in particular, to model the dynamics of genetic 
circadian pacemakers, hormonal cycles and some metabolic pathways. 

In this paper we point out that for Goodwin’s oscillators, where the 
individual reactions have nonlinear (e.g. Mikhaelis-Menten) kinetics, the 
synchronization criterion, obtained by Hamadeh et al., cannot be directly 
applied. This criterion relies on the implicit assumption of the solution 
boundedness, dictated also by the chemical feasibility (the state variables 
stand for the concentrations of chemicals). Furthermore, to test the 
synchronization condition one needs to know an explicit bound for a 
solution, which generally cannot be guaranteed under linear coupling. 
At the same time, we show that these restrictions can be avoided 
for a nonlinear synchronization protocol, where the control inputs are 
“saturated” by a special nonlinear function (belonging to a wide class), 
which guarantees nonnegativity of the solutions and allows to get explicit 
ultimate bounds for them. We prove that oscillators synchronize under 
such a protocol, provided that the couplings are sufficiently strong. 


I. Introduction 

The rhythmicity of many vital processes in living organisms, such 
as the cell division, blood pulse and breathing, diurnal sleep and 
wake cycle, are controlled by genetic and other biochemical “clocks”, 
or pacemakers, that are typically described by nonlinear systems of 
differential equations with stable limit cycles as their solutions. One 
of the first and most influential models of this type, describing genetic 
oscillators [l]-[4], metabolic pathways in a cell [5] and hormonal 
cycles [1], [6], [7], is known (along with its extensions) as Goodwin’s 
oscillator. For 50 years since Goodwin’s seminal paper [8] this model 
has been attracting intensive attention in applied mathematics. 

A challenging problem concerned with biochemical oscillators is to 
study mechanisms of their synchronization via coupling. Experiments 
and extensive simulations (see [2]—[4] and references therein) show 
that the stable 24h-periodic circadian rhythm is not inherent to 
intracellular genetic oscillators (whose natural periods are spread 
from 20h to 28h) but emerges due to the coupling among them, 
which also facilitates the oscillations’ entrainability by the daylight 
and other environmental cues (“zeitgebers”). 
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The Goodwin oscillator is a special case of a cyclic feedback 
system (CFS), consisting of incrementally output passive blocks. An 
important step in understanding the synchronization mechanism for 
such systems has been done in the recent paper [9], establishing 
an elegant criterion for synchronization of identical CFS under 
sufficiently “strong” linear diffusive couplings. An ensemble of CFS 
gets synchronized if the algebraic connectivity of the (weighted) 
digraph, describing the coupling between the systems, exceeds the 
incremental passivity gain of the CFS. A critical observation is that 
this gain depends on the secant gains of the constituent subsystems. 

As will be shown, the criterion from [9] not only adopts an 
implicit assumption of the solution’s boundedness, but in fact requires 
to find the bounds explicitly. In order to apply this criterion, one 
needs to estimate the incremental passivity (or secant) gains of the 
blocks constituting the CFS. Synchronization is guaranteed only when 
these gains are finite, except for that of the leading block, since 
otherwise the minimal coupling strength, required to synchronize 
oscillators, becomes infinite. As will be discussed in Section II-B, 
the chemical reactions with linear kinetics correspond to the blocks 
with finite secant gains. However, nonlinear (e.g. Mikhaelis-Menten) 
kinetics, typically arising in models of enzymatic and other biochem¬ 
ical reactions [10], [11], lead to the infinite passivity gain of the 
correspondent block. This gain becomes finite only for solutions, 
confined to some bounded set, and to estimate the gain, one has 
to find this set or, equivalently, explicit bound for the solution. For a 
linear diffusive protocol establishing such bounds is a non-trivial, and 
in fact open problem. This hinders application of the criterion from 
[9] to Goodwin’s biochemical oscillators with Mikhaelis-Menten 
nonlinearities, modeling e.g. the genetic circadian clocks [2]—[4]. 

In this technical note, we propose a modification of the algorithm 
from [9], combining the usual diffusive coupling with a nonlinear 
“saturating” map, which, similar to the linear protocol from [9], 
guarantees non-negativity of the solutions but, additionally, provides 
an explicit upper bound for the solutions. Under some technical 
assumptions, we prove that the ensemble of CFS’s synchronizes, and 
find explicitly the margin for the coupling strength. Unlike linear 
coupling protocols, the “saturated” protocol also guarantees non¬ 
negative control input which can be important when such an input 
stands for some chemical concentration (e.g., the models of circadian 
oscillators from [2]—[4] treat the input as the concentration of a 
neurotransmitting polypeptide in the extracellular domain). 

The main contribution of the paper is twofold. First, we point out 
some limitations of the synchronization criterion from [9], concerned 
with the necessity to prove the solution boundedness and estimate the 
incremental passivity gains. Second, we develop the approach from 
[9] to address “saturated” protocols, providing synchronization of 
Goodwin-type oscillators with nonlinear reactions’ kinetics. Dealing 
with a more general class of Goodwin’s oscillators, our result 
inevitably inherits two basic limitations of the incremental passivity 
approach [9], [12] and is confined to identical oscillators and diffusive 
couplings (the input of each oscillator depends on the deviation 
between its own and neighbors’ outputs). The results of this paper 
can be applied e.g. to synchronization of synthetic oscillator networks 
(see e.g. [11] and references therein), where individual oscillators and 


2 


synchronization protocols are artificially engineered. 

Note that such oscillator networks as the main circadian pacemaker 
in mammals consist of heterogeneous cells that are coupled non- 
diffusively (being, in fact, an example of pulse-coupled network). 
A simplified continuous-time model for such a network, proposed 
in [2], [3], employs a non-diffusive mean-field coupling. Unlike 
the diffusive protocols, under mean-field coupling the inputs of 
oscillators are identical (depending on the average concentration of 
neurotransmitter, released by individual cells) and do not vanish 
as the oscillators get synchronized. Synchronization of oscillators 
under mean-field couplings and more complicated “nearest-neighbor” 
coupling rules [4] remains a non-trivial mathematical problem. 

II. A CLASS OF CYCLIC FEEDBACK SYSTEMS WITH 
INCREMENTAL PASSIVITY PROPERTIES 

We first briefly recall the central concepts of the incremental 
output strict and feedback passivity (iOSP and iOFP) [9], To sim¬ 
plify matters, we confine ourselves to single input-single output 
systems. We will make use the following notations. Given a vec¬ 
tor x = (xi,... ,x n ) T £ IR", we denote ||a:|| = Vx T x and 
||a:||oo = max; \xi\. We define ljv = (1,1,..., 1) £ IR 1 ' 1 , and use 
the symbol Jjv for the identity N x W-matrix. Given a matrix L, let 
HLIloo = max, fc \Ljk \ stand for its max norm; for z = Lx one has 
then ||z||oo < ||R||oo||a:||oo. Let 1R+ = (0;+ 00 ) and 1R + = [0; 00 ). 


such a system H incrementally output feedback passive with gain 7 
in the set Q x , written iOFPfy -1 ,^); if 7 > 0, we call this property 
incrementally output strict passivity in Q x , written iOSP (7 .CL;). 


B. An iOFP criterion for CFS 

Many cyclic feedback systems (CFS), including Goodwin-type 
oscillators, appear to be iOSP or iOFP, provided that all of their 
sub-systems are iOSP. The results of [9] are concerned with CFS 
whose structures are described by the block diagram in Fig. 1. As 
illustrated, the overall system with the input u ex t £ IR and the output 
y\ £ IR consists of n > 1 nonlinear subsystems H;, governed by 


Vi 


lfii(Xi,Ui) 
Qi{Xi,Ui) ’ 


(4) 


Here Xi £ IR 1- *, u; £ IR, yi £ IR are H;’s state, input and output 
respectively, and pi : IR r * x IR —i IR'* and p; : IR r * x IR. —i IR are 
Lipschitz. The cascaded structure of the system imposes that 



Fig. 1: Block diagram of a cyclic feedback system 


A. Definition of the iOFP property 

Definition 1: A system H, whose dynamics obey 


Ui = yi-i, i = 2,...,n 

U 1 — Uextifi) t/n* (5) 


H : 


x 

V 


p{x,u) 

q{x,u) 


( 1 ) 


where x £ FT 1 , u £ IR, and y £ IR. stand respectively for the state, 
input and output of H, is said to be incrementally output strictly 
passive with passivity (or secant) gain 7 > 0, written iOSP( 7 _1 ), if 
a radially unbounded, positive definite function S : IR' 1 —> IR. exists 
such that for any two solutions to ( 1 ), denoted repectively by x + 
associated with y + , u + and x* associated with y, u ', the increments 
A* = x + — x. Ay = y + — y' and A u = it + — u' satisfy 

Ax) < AuAy - 7 _ 1 |At/| 2 . (2) 

The function S is referred to as the incremental storage function. 
More generally, H is said to be incrementally output feedback passive, 
written iOFP( 7 ” 1 ), if inequality (2) holds for some nonzero 7 £ 
IR U {+ 00 } that has been relaxed from being strictly positive. 

In a degenerate case of static input-out map of H, taking the special 
form u 1 —i y — g(u), inequality ( 2 ) simplifies to the condition 

0 < AuAy — — |Ay| 2 , (3) 

7 

i.e. g is non-decreasing and Lipschitz: |<?(ui) — g(u 2 )\ < y\ui — W 2 |. 

Definition 1 deals with the case when system H defined globally, 
that is, x(t) £ IR r may be arbitrary. The dynamics of biochemical 
systems are naturally defined in the positive orthant, and the iOFP 
property for such systems often can be proved in even more narrow 
domains. We say the inequality (2) is satisfied in a set Q x C IR r , 
if it holds for any two solutions x + associated with y + ,u + and x' 
associated with y**,u* as long as x + (t) and x\t) have the property 
that x + (t), x'(t) £ Q x for any t > 0. Here Q x is not necessarily 
invariant and we are only checking those solutions that stay in Q x 
for all t ; the maps p, q are defined on Q x x FT'" and S in (2) is defined 
at least on the set Q x — Q x = {xi— X 2 £ lR r : xi,X 2 £ Q x }, positive 
definite and radially unbounded (if Q x — Qx is unbounded). We call 


Therefore, the dynamics of the overall CFS can be described by 


Xi 

= pi(xi,u ext — y n ), 

yi = 

Ql (^1 ? Uext 

X2 

= P 2 (x 2 ,yi), 

V2 = 

Q2(x2,yi ), 

x n 

— Pn(x n ,yn— l), 

y-n — 

Qn 5 Un — 


An important result of [9, Theorem 1] states that a CFS, composed 
of iOSP blocks H;, is always iOFP with a gain, satisfying the secant 
condition. Namely, if each H; is iOSP( 7 ~ 1 , Q Xi ) (with 7 ; > 0), then 
CFS ( 6 ) is iOFP(— k,Q x ), where Q x = Q X1 x ... x Q Xn and 

- a 1 / 7r\™ 

k>k = - 1-7273 cos- • (7) 

71 V n) 

Theorem 1 in [9] provides 1 a constructive way to find the incremental 
storage function V(Ax) (positive and radially unbounded), such that 

V(Ax) + a||Ar/|| 2 < fc(Ayi) 2 + Ay^Aiiext (8) 

for any two solutions staying in Q x . Here x = [xj,..., xZ ] T , y = 
[y \,..., y n ] T stand for the joint state and output vectors respectively. 

Notice that k < 00 if and only if all the gains 7 ; are finite, 
except for possibly the gain 71 of the leading block. For the blocks, 
constituting the Goodwin-type oscillator, the gain is usually finite 
only in a bounded domain due to presence of the Mikhaelis-Menten 
or the Hill nonlinearity . 2 Basically, the state-space description (4) of 
the block H;, representing a chemical reaction, is given by 

Xi = -fi(xi)+Ui £ FI, yi = gi(xi) £ FI, (i = 1, 2,..., n) (9) 

where the maps /,. g t are strictly increasing and Lipschitz continuous. 

'Formally, this theorem deals only with the case where Q Xi = ]R r * Vi, but 
its extension to CFS that are not globally iOSP is straightforward. 

2 The Mikhaelis-Menten function is the nonlinear function in the form 
f(x) = f\ \ x./(K ’2 + x), and the Hill function is the nonlinear function 
in the form /( x) = /(Ii '2 + x p ), where K\, K 2 ,p > 0 are constant. 






3 


In the case when g; is linear 3 {g((x) = const > 0), it was proved 
in [9], [13] that the subsystem (9) is iOSPUy” 1 ), where 


and output of the jth CFS in the network, j = 1,..., N. Consider 
the control law, forcing the inputs of the coupled CFS in the form 


7i 


= sup 

Xi 


g'i(xj) 

film)' 


( 10 ) 


This claim remains valid also for nonlinear monotone functions / t , gi 
if inf a- g'i{xi) > 0 [13]. The following lemma extends this result to 
the case where /, is not necessarily monotone and shows that the 
condition inf g' > 0 is critical and cannot be dropped. 

Lemma 1: Let Q Xi C IR be an interval where the functions fi,gi 
are C 1 -smooth. For system (9), the following claims hold: 

1) If infa;, g'i(xi) > 0 and ki = inf*, ++ > — oo, then (9) is 

9i\ x i) 

iOFP(fci,C+). Hence, if fi is non-decreasing on Q Xi . i.e. ki > 0, 
then system (9) is iOSPFy” 1 ^^) with 7 , given in (10). The 
incremental storage function S(Ax) can be chosen quadratic. 

2) Let Q Xi = (0; 00 ), fi be globally bounded and gi(x) ~ a + 
bx~ a as x —¥ + 00 , where a, b £ IR and a > 0. Then (9) can 
never be iOSP no matter what passivity gain is chosen. 

Proof: We prove 1) first. For notational simplicity, we drop the 

subscript i in fi,gi, Xi, xn, ki,Q Xi throughout the proof. By definition 
of k, one has f- (x) — kig'fx) > 0 for any x £ Q Xi . For any pair of 
solutions z + = (x + ,y + ,u + ) and 2 * = (x\y\ u+ x + (f),x+) £ 
Q x , let A 2 = (Ax, Ay, Am) = z + — z*, Af(t) = f(x + (t)) — 
/(*!(()) and A g(tf = g(x + (t)) —g(x\t)). Applying the mean value 
theorem to / — kg, we know that A/(f) — kAg(t) = (/' (9(f)) — 
kg'(9(t)))Ax(t) for which 9(t) £ Q x lies between x + (t) and x+), 
and hence (Af(t) — kAg(t))Ax(t) > 0 and AfAx > kAgAx = 
k.AyAx. Therefore, we have 


AuAx — kAyAx = Ax Ax + AfAx — kAgAx > Ax Ax. 
Taking S(Ax) = ^-|Ax| 2 , where e := inf g(x) > 0, one has 

Am A y — kAy 2 = (Am Ax — kAy Ax) > eAxAx = S(Ax), 

and thus (9) is iOSP(fc,£/ z ), which proves statement 1). 

We prove 2) by contradiction. Suppose on the contrary that 
(9) is iOSP(fc, Q x ) with the storage function S(Ax). Let M = 
sup | f(x) | < 00 , and x, x + be a pair of solutions under the inputs 

X 

v) (t) = M + (7 + l)f 7 and M + (f) = 2M + e + u+) respectively, 
where (7 + l)a > 1 and e > 0. Consequently, x' > (7 + l)t 7 so 
that x*(t) > f 7+1 + x^(0) as t —» 00 and A x(t) > e. Therefore, 
Ay(t) is in the order of as t —> 00 , and thus AuAy and 

\Ay\ 2 are summable functions. By integrating (2), one obtains that 

r T 

S(Ax(T)) - S{ Ax(0)) < / (AuAy - k\Ay\ 2 )dt < 

Jo 

POO 

< / (|Am|| Ay\ + |fc|| Ay\ 2 )dt < 00 for all T > 0, 

Jo 

which contradicts the fact that S(Ax(T)) —> 00 as T —> 00 . ■ 


III. Synchronization of diffusively coupled CFS 

The iOFP property of the CFS (6) allows to prove synchronization 
in a network of N > 1 identical CFS, where the couplings are 
described by a weighted, strongly connected, balanced graph with 
the Laplacian matrix L £ lR JVxJV . Let u 3 ext , x 3 = [x-j T ,... ,x 3 n T ] T 
and y 3 = [y{,... ,yi] T denote respectively the external input, state 

3 The linearity follows implicitly from the assumption that the integral in 
[9, eq.(8)] is well-defined. 


u ext (t) = —LYi(t), 


(ID 


where U ext = [u \ xt , ... , u* xt ] T and IT = [y\, y?} T . 

The result of [9, Theorem 2] shows that protocol (11) synchronizes 
the outputs of N CFS’s if the coupling is sufficiently “strong”; its 
“strength” is bounded below by the algebraic connectivity of the 
graph if L — L T or, generally, by the second smallest eigenvalue 
A 2 of the matrix L+ Jf . Precisely, let each CFS satisfy (8) for some 
a > 0 and k > 0 and be limit set detectable [9], If A 2 > k, then 
any bounded solution of the coupled CFS achieves synchronization 

lim \x{(t)-x k i(t)\=0 Mi = 1,..., n; Mj, k = 1,...., N. (12) 

£—>• + 00 


Remark 1: Theorem 2 in [9] claims a more general result stating 
that synchronization is achieved without the assumption on bounded 
solutions; however, as discussed below, this assumption is implicitly 
required in its proof when appealing to the LaSalle invariance 
principle. For special types of oscillators, e.g. the Lur’e system 
with sector nonlinearity, the solution’s boundedness is ensured by 
the input-to-state stability property of the individual system [12]. 
However, in general the technique to drop it remains elusive if not 
impossible. 

The following theorem extends Theorem 2 in [9] to the case where 
the iOFP property holds only in some domain; its proof, following 
the line of the proof from [9], demonstrates, in particular, that the 
boundedness assumption is essential. 

Theorem 1: Consider a system of N > 1 identical limit-set 
detectable [9] CFS (6), satisfying (8) with some a, k > 0 in some 
closed domain Q x C IR ri+ "' +rn . Suppose the CFS are coupled 
together through the protocol (11) with A 2 > k. Then any bounded 
solution of the closed-loop system, such that x 3 ( t ) £ Q x Mt> to for 
j = 1,..., N and some to > 0, asymptotically synchronizes (12). 

Proof: As before we use x 3 £ lR riH l_r " to denote the state 
of the jth CFS, and now let £ = [x lT ,..., x nT ] t and £ = 
[Y lT ,... , Y nT ] t be the state and output of the overall networked 
system respectively. For the incremental storage function V of each 
individual CFS, which satisfies (8), and for 1 < j,m < N, let 
Vj,m(i) = V(x m -x 3 ) and S+) = A. J2* m=1 Vj,m(£). Substitut¬ 
ing solutions (x p ,u p ,y p ) and (x q , u q , y q ), where x p (t),x q (t) £ Q x 
for t > to, into (8), the following condition is valid as t > to: 

V(x p - x q ) +a\\y p -y q || 2 < k(y p -y q ) 2 + (y p - y q )(u p xt - u q ext ). 


By summing up these inequalities over all p, q and introducing the 
projector n = In — -^IjvIjv, one arrives at [9] the following 

s(t) < -a II(n ® LOCH 2 + (fclinTTf + (m\) T nu RXt ). (13) 

Using (11), one easily finds that UU ex t = —HLY\ = — UIYi and 
hence (niT) T n[/ ex t < —A 2 ||HVi|| 2 ; therefore (13) implies that 

S(0 < -a||(n® 7 n )C(f)|| 2 < 0 Mt > to- (14) 

Let £(t) be a bounded solution with x 3 (t) £ Q x and M = {£ = 
(x 1 ,...,x JV ) £ Q x : S(0 < S(£(to))}. Due to (13) one has 
£(t) £ M and hence the closed set M contains the cu-limit set 
of £(•). Thanks to the LaSalle invariance principle, the solution £(f) 
converges to the maximal subset of M, where S = 0 and hence 
y 1 = ■ ■ • = y N due to (14). The limit-set detectability assumption 
entails now synchronization of the state vectors. ■ 

Note that a widely used version of LaSalle invariance principle [14] 
requires the Lyapunov function S, along with S, to be defined on a 
compact invariant set A4, and guarantees that any solution starting 
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in M converges to the maximal set where 5 = 0. However, original 
versions of LaSalle’s invariance principle [15], [16] are applicable to 
any bounded solution and guarantee that S = 0 on its cc-limit set, 
provided that S and S are well defined in the vicinity. Assumption 
of compactness and invariance of M automatically entail the latter 
condition, as well as boundedness of any solution starting at M. 
Without this assumption, LaSalle's invariance principle can still be 
applied, but the extra condition of boundedness is then unavoidable. 

To prove synchronization of CFS under linear balanced protocol 

(11) , using Theorem 2 in [9] or more general Theorem 1, one has first 
to establish the iOFP property in some domain Q x . As follows front 
Lemma 1, the relevant passivity gains 7 i of the subsystems H, can be 
infinite or even undefined, unless the corresponding state variables Xi 
are confined to some bounded domains Q Xi . So the restriction x(t) £ 
G x , imposed to apply the iOFP property ( 8 ). requires to find some 
explicit bound for the solution. Even if Q x can be unbounded (like in 
the example from [9]), the criterion still guarantees synchronization 
only for bounded solutions. Using (14), deviations x p — x q are shown 
to be bounded, entailing boundedness of the states x 3 ( t ) under input- 
to-state stability assumptions (which hold e.g. for Lur’e-type systems 

[12] ). However, proving the solution boundedness for general CFS, 
coupled via a linear protocol ( 11 ), remains a non-trivial problem. 

To cope with this problem, we replace the linear protocol (11) 
with a nonlinear one, providing sufficiently small and non-negative 
control inputs u 3 ext , that is, 0 < u 3 ext (t) < Mo, where Mo is some 
known constant. Under such a constraint, one often can localize the 
solution in a domain where the incremental passivity gains of all 
the subsystems are known and finite. Relevant sufficient conditions 
for this, dealing with Goodwin's oscillators, will be discussed in 
Section IV. In fact, the input restrictions are often dictated by the 
biological feasibility, e.g. in some models of coupled circadian clocks 
[2]—[4] the oscillators’ inputs stand for the concentrations of the 
neurotransmitter in extracellular media. In this technical note, we 
do not aim to examine the model from [2]—[4] itself, which consid¬ 
ers mean field couplings. Instead, we propose a diffusive coupling 
protocol similar to ( 11 ), but employing a non-negative “saturating” 
nonlinearity, which guarantees the input constraint and thus entails the 
solution’s boundedness. Meanwhile, the protocol constructed below 
provides synchronization under sufficiently “strong” coupling, and 
the minimal sufficient strength may also be explicitly estimated. 

The algorithm we propose is as follows 

u 3 ext {t) = g 0 (cv 3 (t)), V(t) = = -LYi(t). (15) 


A constant c stands for the coupling gain; here, to emphasize the 
effect of the coupling strength, we have intentionally added c that 
has been implicitly incorporated into the entries of L in (11) as is 
done in [9]. The function go : IR —» [0; + 00 ) is bounded, saturating 
the inputs at a prescribed constant Mo = sup„ gH go (u). The auxiliary 
inputs v 3 ( t ) £ IR are introduced to emphasize the similarity between 
the protocols (11) and (15): in fact, the system of N CFS’s (6), 
coupled through the protocol (15), may be considered as a collective 
of appropriately modified CFS’s, coupled linearly. 

Hereinafter, we assume the following assumption to be valid. 
Assumption 1: The function go(-) is smooth, globally bounded and 
strictly increasing, hence g'o(v) > 0, Vv £ IR. Additionally, u(s) = 
inf g'o{v) decreases at infinity more slowly than linear functions, 

\v\<s 

i.e. v(s) —» 0 yet |s|^(s) —» +00 as s — t ± 00 . 

Assumption 1 is satisfied by a wide class of functions, e.g. 


ffoO) = 


Mo 


1 + 


[upsignu 

T+TuF 


with Mo > 0, p £ (0,1). (16) 


Note that Mo = sup</o(v) may be chosen as small as possible. 

u€K 



Fig. 2: Auxiliary cyclic feedback system with saturated input 


The next result shows that the modified protocol (15) synchronizes 
the systems (6), provided that the coupling is sufficiently strong and 
the solution stays in some compact set; the crucial difference with the 
linear protocol (11) is that the existence of such a set attracting the 
solutions may often be proved by choosing </o(') sufficiently small. 

Theorem 2: Suppose that Assumption 1, and the assumptions of 
Theorem 1 hold, where Q x C IR ri + -+ r '“ is a compact set and 
p\(x\,U\) = pi(*i). Then for sufficiently large gain c, any solution 
of the closed-loop system such that x p (t) £ Q x Vt > to for all p = 
1,..., N and some to > 0 gets synchronized (12). Synchronization 
is implied by the following inequality, which holds as c —¥ +00 

ciz(c||I/||t/*)A 2 > k, y» = max{|pi(a;i)| : x £ Q x }- (17) 

Proof: Along with the original CFS (6), consider a modified 
system with a new input v(t), which obeys (6) and the addi¬ 
tional equation u ex t = cgo{cv) as shown in Fig. 2. The network 
of CFS (6), coupled via protocol (15), is now equivalent to the 
network of N “augmented” systems, coupled via (11). One may 
easily notice that if x p (t) £ Q x then | 2 / p (t)| < y * and hence 
|u p (t)| < ||F'||oo)/* due to (11); this implies that 0 < u p xt (t ) < 
m* = max go- Due to mean value theorem, for any two solutions 
one has A u ex t(t) = cg' 0 {c9(f))Av, where |#(t)| < ||L||t/*, and thus 
0 < m = ciz(c||Z/||oo?/*) < cg'o(c9) < M = cmaxpo(u) and thus 

d€IR 

the right-hand side of (8) is not greater than M(m~ 1 k\Ayi(t)\ 2 + 
Av{t) Aj/i(f)). This ensures that the augmented CFS also satisfies 
(8), where u ex t, a and k are to be replaced with respectively v, 
a = AT~ 1 a and k = k/m. Synchronization now follows 4 from 
Theorem 1 since A 2 > k due to (17). ■ 

IV. Applications in coupled biochemical oscillators 

In this section, we discuss how Theorem 2 allows to estimate the 
coupling strength, needed to synchronize biochemical oscillators of 
the Goodwin type, governed by the equations 


Xl 

= -fl(xi) + (Uext - Vn) 

i, yi = gi(xi), 

X2 

= -/2(x 2 ) + yi, 

2/2 = 52 (x 2 ), 

Xn-1 

= fn — l (Xn—1 ) ~\~ yn—2i 

2/n— 1 = Qn.— l{X: 

Vn 

= 9n(yn — l)- 



Therefore, the Goodwin-type oscillator is an example of the system 
(6) shown in Fig. 1, where the blocks Hi,..., H n _i obey equations 
(9) with Mi = Uext — y-n, U 2 = yi, ■■■, u n = 2/n-i, and the block 
H„ is static: y n {t) = g n (u n {f)). The classical Goodwin's model 
[8] corresponds to the case where fi are linear and g n is the Hill 
nonlinearity. We emphasize that the system operates in the positive 
orthant, namely Xi > 0, 1 < i < n — 1. 

Consider now a network of N identical oscillators (18), coupled 
via the distributed protocol (15), where go satisfies Assumption 1. 

Denoting Mo = max go (v), we adopt the following assumption. 

4 Formally, Theorem 1 was formulated for systems (6). However, its proof 
employs only the inequality (8) and obviously remains valid in spite of the 
additional static block Hq. 
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Assumption 2: The functions fi,gi : IR+ 1R+, i = 1 ,,n — 
1. go : IR H > [0,M o \,g n : IR+ [— M n ,0] are smooth and strictly 
increasing (/,■. g' > 0). Additionally, the maps f i satisfy the condition 

$;(0) = -oo, T>i(+oo) = +oo, $;(*) = [ -jy-t,x>0. (19) 

J i /i(s) 

Condition (19) holds, for instance, for linear functions fi(x) = 
aiX, a,i > 0 and Mikhaelis-Menten type nonlinear functions. Under 
Assumption 2, any solution of the cyclic feedback system with 
saturated input, starting strictly inside the positive orthant, remains 
positive and, under additional assumptions, is ultimately bounded. 

Lemma 2: Let Assumption 2 hold and u ex t(t) £ [0; Mo] is defined 
for t 6 A = [0; /?). Then for any initial condition Xi( 0) £ IR+ (i = 
1 ,..., n— 1 ) the solution xi(t ),..., x n -i(t), yi(t),..., y n (t) exists 
on A and remains positive Xi(t) £ IR+Vf £ A. If (3 = oo and the 
functions hi. given by the recursion hi = /j -1 , h .2 = ff 1 og 1 ohi ,..., 
h „-i = f~l i °g n -2 °hn- 2 , are well defined on [0; M + £o), where 
M = Mo + M n and £o > 0, then the solution is ultimately bounded 

lim Xi(t) < Xi Vi = 1,..., n — 1 Vxj(O), (20) 

t—y oo 

where Xi = hi(M) are independent of the initial conditions. 

Proof: We now prove that x(t) exists and is positive on A. 
Since Xi(0 ) > 0, because of continuity, there exists a maximal 
interval A' = [0; f3') C A, such that Xi(t ) > 0 for all t £ A' 
and i = 1,... ,n — 1. From Assumption 2, for t £ A' it always 
holds that Ui(t) > 0 since y n (t) < 0, tti = u ex t — Vn, U 2 = yi, 

Un = y-n—i. Therefore 4-T>i(*i(t)) = Xi{t) / fi(xi(t)) > —1, which 
implies that Xi(t) > $£ ('Fi(O) — f) > 0 for any t £ A'. Hence, the 
solution cannot escape from IR+ within A\ Furthermore, if (3 1 < oo 
then Xi(t) > Si > 0 and hence fi{xi(t )) > Vi > 0 as t £ A'. This 
implies that ( 07 (f)) and hence Xi(t) are bounded from above, i.e. 
the solution cannot grow unbounded in finite time. This, according 
to the definition of A', implies that A' = A. 

To prove (20), we show first that for any solution Xi(t) > 
0 ,Ui(t) > 0 of the subsystem (9) the following implication holds: 

Pi = lim m(t) < sup fi(x) =7 lim 07 (f) < ff 1 (pi). ( 21 ) 

t —F + oo x>0 t—>oo 

Indeed, let 8, e > 0 be so small that pi + 8 + £ < sup From the 
definition of the upper limit, a number To exists such that Ui(t) < 
Pi + 8 for all t > To. From the facts that 07 (f) < —e if f > To and 
Xi(t ) > £ = ff 1 (pi + 8 + s ) fi(xi(t)) > pi + 8 + s > Ui(t) + s, 
we know the following two statements hold: (i) if Xi(Ti) < ^ for 


non-trivial problem, which is beyond the scope of this technical note. 
However, the following simple lemma shows that in practice (20) 
holds as u ex t(t) is sufficiently small, provided that the oscillators 
have a globally stable attractor (e.g. limit cycle). 

Lemma 3: Suppose that any solution of the system (18) with 
Uext = 0, starting at a compact K C IR" -1 , converges to some 
attractor Kq (2 IntK: dist(x(t), A'o) — > 0 as t. —» oo. Then for any 
8 > 0 there exist £o = £o(A'o, A') > 0 such that any solution of (18), 
starting at *(0) £ A' and associated with input 0 < u ex t(t ) < £o, 
converges to the attractor's 5-neighborhood. Precisely, there exists 
To = To(8, K, A'o) such that dist(x(t), Ko) < 8 as t > To. 

Proof: Without loss of generality, let 8 > 0 be so small that 
Ks = {a: : dist(x(t), A'o) < 5} C K. Since K is compact, there 
exists To such that x(t) £ Ks as t > To under u ex t = 0 for any 
solution, starting at *(0) £ K. Hence, if 0 < u ex t(t) < £oVf £ 
[0; 2To] and £o > 0 is sufficiently small, one can guarantee that 
x(t) £ Ks C K at least for t £ [To;2To] independent of the initial 
condition in K and concrete u ex t. Applying this for *(0) = x(To) £ 
K and shifted input u ex t{s) = u ex t(To + s), one shows that x(t) £ 
Ks C K as t £ [2To; 3To], and so on. ■ 

Assumptions of Lemma 3 in general hold for biologically realistic 
models, where the oscillating concentrations of the reagents are 
confined to some (roughly known) intervals and the limit cycles are 
found experimentally or via numerical simulations. Lemma 3 states 
that knowledge of the attractor allows to estimate the solutions of the 
network of coupled oscillators, provided that the control inputs are 
sufficiently small. However, unlike Lemma 2, Lemma 3 does not give 
the explicit dependence between x l and the value of ATo max n e a;t, 
but only allows to find the limit of x l as Mo —> 0 (and x(0) £ A'). 

We now return to the dynamics of coupled CFS (18) under 
“saturated” protocol (15) and prove that if the ultimate boundedness 
(20) holds under “weak” non-negative inputs u ex t{t) (for instance, 
condition from Lemma 2 or Lemma 3 is valid), then (15) guarantees 
synchronization for small Mo is small, if the coupling is sufficiently 
strong: c > c* = c*(xi) (and c* can be found explicitly). 

For convenience, we introduce the following assumption. 

Assumption 3: For some set Q C R" -1 there exist such Mo > 0 
such that under any input u ex t(t ) £ [0; Mo\ Vf > 0 solutions of (18), 
starting at *(0) £ Q, satisfy (20), where the bounds x 1 are uniform 
over all *(0) and u ex t{-). 

Assumption 3 can be provided, for instance, by the conditions 
from Lemma 2 or Lemma 3. If it holds, any solution which starts 
at Q enters in finite time the hypercube B e = {x : 0 < Xi < xf = 


some Ti > To, then Xi(t) < £ Vf > Ti, and (ii) such a Ti necessarily 
exists; that is, if Xi(To) < £, one can take Ti = To, otherwise, 
Ti < ( Xi(To ) — £)/e is the first time instant after To at which 
Xi (Ti) = £. Hence, lim Xi(t ) < f~ 1 (pi + 8 + e), front which 

t—y “I - OO 

( 21 ) follows by passing to the limit 8, e —» 0 . 

Then (20) can be easily proved using (21). From the assumptions, 
the input of the block Hi is given by ui(t) = u ex t(t) — y n (t) < M. 
Therefore, lim xi(t) < ff 1 (M) = hi(M) and hence lim yi(t) < 

t—too t—too 

gi o hi{M). Invoking (21) for the second block H 2 with input U 2 = 
2 / 1 , we obtain lim X 2 (t) < ff 1 ° gi o hi(M) = h 2 {M). Iterating 

t—t OO 

this procedure for H 3 ,..., B„_i, the inequalities ( 20 ) are proved. ■ 
Lemma 2 gives only the simplest condition of “restricted” input- 
to-state stability (ISS) [17], that is, the existence of explicit ultimate 
bounds for the state vector of CFS provided that its input is suffi¬ 
ciently small and positive. This condition appears to be conservative 
for some Goodwin-type oscillators, as will be discussed below. It 


Xi + £, 1 < i < n} and remains there. Since the blocks Hi are 
iOSP(( 7 i ) -1 , [0; Xi]) for i = 1,... ,n by virtue of Lemma 1, where 
the secant gains are given respectively by 


e A 

7 i = max —7—7, 
f'(Xi) 


i <n, 7„ 


max 

y£[0,g n -l(.x 


As follows from Theorem 1 in [9] (see discussion in Section II-B), 
the CFS (18) is iOFP(-fc £ ,H e ) and, moreover, the inequality holds 


14(Aa:(f))-l-a s ||Aj/(f)|| 2 < k e (Ayi(t)) 2 + Ayi(t)Au ex t(t), (23) 


for any two solutions of (18), staying in B e . Here k e = — ^ + 
7I73 ... 7^. Notice that £ > 0 can be arbitrarily small, independent 
of the initial condition. If x{t) £ B e , then 0 < t/i(f) < y% = gi{x\). 

Theorem 3: Let N identical CFS (18) satisfy Assumptions 2 and 
3 and be coupled via (15). Suppose that go satisfies Assumption 1, 
go(v) < MqVv £ R and the coupling is sufficiently strong 


can be further refined (with tightening the bounds x l ) by using the 1 0 0 0 

monotonicity-type arguments from [18], [19]. To establish the ISS a/ ( c ll^ / ll?A )^2 > ~^o+7273 •••7m 
property for general CFS (with explicit bounds x') remains an open 


0 A e . 0 A s , 

7 i = it \e=o,y» = y,U=o- 


( 24 ) 
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Under such coupling, any solutions (x° (t))j = i, starting at a: p (0) £ Q, 
remain in the positive orthant xl(t) > 0 and synchronizes ( 12 ). 

Proof: Since 7 ? is continuous at e — 0, inequality (24) remains 
valid after replacing 7 ° 1 —>■ 7 ? and y° i-> yf when e > 0 is small. 
The claim now follows from (23) and Theorem 2. ■ 

V. Example: synchronization of circadian clocks 

In this section we demonstrate our synchronization criterion for 
oscillators, describing the cellular circadian clocks [2], [3]. The main 
circadian pacemaker in mammal is controlled by the neurons of SCN 
(suprachiasmatic nucleus, a zone in hypothalamus). Within each cell 
(indexed 1 through TV), a clock gene mRNA (X) produces a clock 
protein (Y) which, in turn, activates a transcriptional inhibitor (Z), 
closing a negative feedback loop [3]; their dynamics are given by 

K n X 

* ^ = Vl K? + (Z*) n ~ U2 1< 2 + + Uext 

Y\t)»ks x '-u< J £^ (25) 

Z\t) = k 5 Y i -v 6 * = 1, 2,..., TV 

Here u' ext are some external inputs. In [2]—[4] the networks with 
mean-field couplings are considered, where u\ xt = ... = u™ xt is 
a common input, which is positive, bounded and depend on the 
average concentration of neurotransmitting peptide in the extracel¬ 
lular domain, depending in its turn on X 1 , ..., X N . We consider 
synchronization of oscillators (25) under diffusive protocol (15). 
Unlike the mean-field control, which remains oscillatory when the 
synchronization is established, under protocol (15) the inputs stabilize 
at the constant value ul xt (t ) —> const = (?(0)V* as t —> 00 . 

It is confirmed experimentally [2] that an individual circadian 
clock has a stable limit cycle in the positive octant X, Y, Z > 0, 
and the corresponding oscillation period lies between 20 and 27 
hours. This means that for realistic sets of parameters in (25) 
Lemma 3 and Theorem 3 work, stating that for any compact set 
K C IR+ one can find Mo = Mo(K) > 0 (sufficiently small) and 
c > 0 (sufficiently large), such that the protocol (15) synchronizes 
oscillators (25) starting at (X*(0), Y l (0), Z l ( 0)) E K. For a special 
set of parameters n, hi, I\i, Ui, found in [3], Lemma 2 is applicable 
which gives explicit estimates for Mq and explicit bounds for the 
solutions. We simulated dynamics of a more complicated model [2], 
where Lemma 2 is unapplicable (in the notation of Lemma 2, the 
function hi = ff 1 is not defined on [0; M n ]) and Mo, c are to be 
found numerically. Note that the incremental passivity of (25) in the 
whole positive orthant does not follow from the criterion in [9], unlike 
the Goodwin oscillator in the example from [9], so synchronization 
of CFS (25) under strong linear couplings remains an open problem. 

We simulate the dynamics of TV = 10 all-to-all coupled oscillators 
(25) with the parameters from [2]: v\ = 0.7nM/h; AT = 1«M; 
n = 4; U2 = 0.35 nM/h; AT = nM; k 3 = 0.7/h; U4 = 0.35nM/h; 
AT = 1 nM; fcs = 0.7 /h; vq = 0.35 nM/h', AT = InM , which 
correspond to the oscillation period ~ 23.5 h. We choose go(v) in 
the form (16), where Mo = 0.0005 and p = 0.9. We simulate the 
dynamics for c = 0, c = 1, c = 10 and c = 100. Oscillators are not 
synchronous for c being small, however, the synchronization emerges 
as c increases, confirming thus Theorem 3. 

VI. Concluding remarks 

Theorem 3 of this note has shown that a similar conclusion in 
comparison to that of [9] holds even when the coupled biochemical 
oscillators are under input constraints, imposed by the requirements of 
biological feasibility and necessity to explicitly estimate the coupling 



Time (h) Time (h) 



c = 10 c = 100 

Fig. 3: Dynamics of the mRNA levels (X 1 ) for c = 0, 1,10,100. 


gains. To satisfy these constraints, we combine the linear coupling 
protocol from [9] with a “saturating” nonlinear block. We have proved 
that strong diffusive couplings can get coupled CFS-type oscillators 
synchronized when the saturation nonlinearity of the oscillators’ 
inputs belongs to the identified class. Our proof is based on the 
synchronization criterion from [9], extended to the systems with 
additional saturated block. The result may be extended to the CFS 
coupled through outputs yk, as considered in [20], The techniques of 
quadratic constraints, used in our recent paper [21], allow to extend 
our results to some other types of “saturated” protocols, where not 
only control inputs, but also outputs (or their deviations) are saturated. 

The results of our paper can be applied e.g. to networks of 
synthetic biochemical oscillators [11] where the couplings between 
the individual oscillators are artificially engineered. However, as 
has been reported by biochemists and biophysicists, the couplings 
between many natural biochemical oscillators, in particular neurons 
of the circadian pacemakers, are in general not diffusive [2], Hence, 
we are studying models for biochemical oscillators under mean field 
coupling [2], [3] or more complicated nearest-neighbor couplings [4], 
regulated by the concentrations of neurotransmitting polypeptides. 
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